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ABSTRACT 

We report the discovery of faint very high energy (VHE, E > 100 GeV) y-ray emission from the radio galaxy 
Centaurus A in observations performed with the H.E.S.S. experiment, an imaging atmospheric Cherenkov 
telescope array consisting of four telescopes located in Namibia. Centaurus A has been observed for more than 
120 h. A signal with a statistical significance of 5.0 cr is detected from the region including the radio core and 
the inner kpc jets. The integral flux above an energy threshold of ~250 GeV is measured to be 0.8 % of the flux 
of the Crab Nebula (apparent luminosity: L(>250GeV)» 2.6 x 10 39 erg s _1 , adopting a distance of 3.8 Mpc). 
The spectrum can be described by a power law with a photon index of 2.7 ± 0.5 sta t ± 0.2 sys . No significant 
flux variability is detected in the data set. However, the low flux only allows detection of variability on the 
timescale of days to flux increments above a factor of ~ 15 - 20 (3 cr and 4<x, respectively). The discovery 
of VHE y-ray emission from Centaurus A reveals particle acceleration in the source to > TeV energies and, 
together with M 87, establishes radio galaxies as a class of VHE emitters. 
Subject headings: galaxies: active — galaxies: individual (Cen A) — gamma rays: observations 
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1. INTRODUCTION 

Centauru s A (Cen A) is the nearest active galaxy (for a re- 
view see Israel| 1998] >. It is classified as a FR I radio galaxy; 
these are thought to be the parent population of BL Lac ob- 
jects ( |Urry & Padovanip9 95 ). At radio wavelengths rich jet 
structures are visible, extending from the core and the inner 
pc and kpc jet to giant outer lobes with an angular extension 
of 8° x 4°. The inner kpc jet has also been detected in X- 
rays, reve aling a complex structure of bright knots and diffuse 
emission ( Kraft et aL|2002| l. The angle of the jet axis to the 
line of sight is estimated to be 15-80° (see e.g. |Horiuchi et al. 
|2006| and refe rences therein). With a distance of 3.8 Mpc 
(Rejk uba|2004| 1' corresponds to ~ 1.1 kpc projected length) 
even the inner jet structures are resolvable with the angu- 
lar resolution of current experiments in the very high energy 
(VHE;£ > 100 GeV) regime 34 . The elliptical host NGC 5128 
features a dark lane, a thin edge-on disk of dust and young 
stars, believed to be the remnant of a merger. Recent esti- 
mates for the mass o f the central superma ssive black hole give 
(5.5 ± 3.0) x 10 7 M Q ( |Cappellari et al.|2008| . The kpc-scale jet 
and the active nucleus are confirmed sources of strong non- 
thermal emission. In addition, more than 200 X-ray point 
sources with an integrated luminosity of Lx > 10 38 ergs _I 
are established to be associated with the host galaxy ( |Kraft| 
et al.|200l| . Recently, [Croston et al] ( |2009| l reported the de- 



H.E.S.S 



tection of non-thermal X-ray emission from the shock of the 
southwest inner radio lobe from deep Chandra observations. 

Cen A was detected at MeV to GeV energies by all in- 
struments on board the Compton Gamma-Ray Observatory 
(CGRO) in the period 1991 - 1995, revealing a peak in the 
spectral energy distribution (SED) in vF v representation at 
~ 0.1 MeV with a maximum flux of about ~ 10~ 9 erg cirT 2 s _1 



(Stein le et al.|1998) . |Steinle et aTT] ( |1998| ) also reported vari- 
ability of the Cen A source, especially pronounced at 10 MeV, 
while Sreekumar et al. ( 1999) found the EGRET flux was sta- 
ble during the whole period of CGRO observations. A tenta- 
tive detection of Cen A (4.5<x) at VHE during a giant X-ray 
outburst in the 1970's was reported by Grindlay et al. ( 1975 ). 
Subsequent VHE observations made with the Mark III ( |Car-| 
Iraminana et al.|1990|l, JANZOS (|Allen et al.|1993| l, CANGA- 
ROO dRowell et al.|| 19991 |Kabuki et al.||2007| l, and H.E.S.S. 
(Aharonian et a l.|2005| l experiments resulted in upper limits. 

Cen A has been proposed as a possible source of u ltra-high 
energy cosmic rays (U HECR; E > 6 x 10 19 " 20 eV) ( |Romero 



et al.| 1996| but see also |Lemoine| 2008). Recently, the Pierre 



Auger Collaboration reported an anisotropy in the arrival di 
rection of UHECR ( |Abraham et "alj|20071 12008] >. While a 
possible correlation with nearby active galactic nuclei is still 
under discussion, it has been pointed out that several of the 
events can possibly be associated with Cen A (e.g. Gorbunov 
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Fig. 1. — Smoothed excess sky map centered on the Cen A radio core 
(cross). Overlaid contours correspond to statistical significances of 3, 4, and 
5<x, respectively. The inlay in the lower left corner shows the excess expected 
from a point source (derived from Monte Carlo simulations). The integration 
radius is 0.1225° and the map has been smoothed with a two-dimensional 
Gaussian of radius 0.02° to reduce the effect of statistical fluctuations. The 
cosmic-ray background in each bin is estimated using the template back- 
ground method iRowell 2003 1. 

|et al.|2008||Moskalenko et al.|2008||Kacherriess et al.|2008] >. 

Until now, the only firmly established extragalactic VHE y- 
ray source with only weakly be amed emission is the giant ra- 
dio galaxy M 87 ( [Aharonian et al.|2003[|2006a[ l. M 87 showed 
strong flux outbursts in the VHE regime with timescales 
on the order of days (Aharonian et al. |2006a Albert et al. 
2008 ), pointing to a characteristic size of the emission region 
< 5 x 10 15 (5cm, corresponding to « 55 Schwarzschild radii 
(M B h = 3.2 X 10 9 M Q , 5: relativistic Doppler factor). Re- 



cently, |Aliu et al.| ([2009) and |Acciari et al.| ( |2009| ) reported 
VHE emission from the direction of the blazar 3C 66 A and 
the radio galaxy 3C66B (angular separation 6')- While |AlIu" 



et al. (2009 ) favor 3C 66 B as the origin of the VHE emission 
in their dataset, |Acciari et al.| ( |2009| l exclude 3C 66 B as the 
origin of their detected emission with a significance of 4.3cr. 

In this Letter the discovery of VHE emission from Cen A 
with the H.E.S.S. experiment is reported, and properties of the 
detected emission and their implications are discussed. 

2. H.E.S.S. OBSERVATIONS AND RESULTS 

The H.E.S.S. (High Energy Stereoscopic System) collabo- 
ration operates an array of four large imaging Cherenkov tele- 
scopes (IACT) for the detection of VHE y-rays, located in the 



Southern Hemisphere in Namibia (Aharonian et al. [2006b ). 
The H.E.S.S. observations of Cen A were performed between 
April 2004 and July 2008. A dead time corrected total live 
time of 1 15.0 h of good-quality data was recorded. The zenith 
angles of the observations range from 20° to 60° with a mean 
zenith angle of ~ 24° . The data were recorded with pointing 
offsets between 0.5° to 0.7° relative to the radio core position, 
to enable a simultaneous estimation of the background us- 
ing events from the same field of view (reflected background) 
( ] Aharonian et al.|2006b[ ). The data were analyzed with a stan- 
dard Hillas-type analysis (Aharonian et al. 2006b i with an 
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Fig. 2. — Optical ima ge of Cen A (UK 48-i nch Schmidt) overlaid with radio 
contours (black, VLA, Condon et al. 1996 1, VHE best fit position with 1 cr 
statistical errors (blue cross), and VHE extension upper limit (white dashed 
circle, 95% confidence level). 
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Fig. 3. — Differential energy spectrum of Cen A as measured by H.E.S.S. 
The line is the best fit of a power law dN/dE = <Dq ■ E~ r to the data (r = 
2.7 ± 0.5 sta i ± 0.2 sys ). The lower panel shows the residual &x~ of the fit. 

analysis energy threshold of ~ 250 GeV for a zenith angle 
of 20°. 

Figure [T] shows the smoothed excess sky map of VHE y- 
rays as measured with H.E.S.S. centered on the Cen A ra- 
dio core position. A clear excess at the position of Cen A 
is visible. A point source analysis, using standard cuts as 
described in |Aharonian et al.| ([2006b ), was performed on 
the radio core position of Cen A, resulting in the detection 
of an excess with a statistical significance of 5.0 cr (calcu- 
lated following Eq. 17 of Li & Ma|1983 330 excess events, 
N 0N = 4199, N OPF = 42513, a = 091). A fit of the 
instrumental point spread function 35 to the uncorrelated sky 
map results in a good fit (chance probability ~ 0.7) with a 
best fit position of 072000 = 13 h 25 m 26.4 s + 4.6* tat + 2.0* yst , 

35 Derived from Monte Carlo simulations. 
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Fig. 4. — Spectral energy distribution of Cen A. Shown are the VHE spec- 
trum as measured by H.E.S.S. (red filled circles), previous u pper limits and 
tentative detect ions in the VHE r egime (purple markers; Grin dlay et al.|1975l 
open diamond; Carraminana et al. 1990 open cross; Allen et al. 1993 filled 
circle; Rowell et al. 1999 open triangle; Aharonian et al. 2005 open circle; 
Kabuki et al.|2U0/| fi lled squares), EGRET measurements in the GeV regime 
(Sreekumar et al.|T9~9 9 grey bow tie), and data from the NASA Extragalactic 
Database (NED) (grey filled circles). 



372000 = -43°0.7' + l.ls tat + 30^, well compatible with the 
radio core and the inner kpc jet region (Fig. [2}. Assuming a 
Gaussian surface-brightness profile, we derive an upper limit 
of 0.2° on the extension (95% confidence level). 

The differential photon spectrum of the source is shown 
in Fig. [3] 36 A fit of a power-law function dN/dE - €>o • 
(£yiTeV)~ r to the data is a statistically good description 
(yVd.o.f. = 2.76/4) with normalization d> = (2.45 + 



1-lstat: 



30:: 



0.52 s 



0.49 sys ) x 10 



13 cm" 2 s- 1 TeV" 1 



and photon index 



r = 2.73 + 0.45 sta t + 0.2 sys . The integral flux above 250 GeV, 

: (1.56 ± 
0.8% of 



taken from the spectral fit, is <t>(£ > 250 GeV) 
0.67 stat ) x 10~ 12 cm -2 s _1 , which corresponds to 



the flux of the Crab Nebula above the same threshold ( |Aha-| 
ronian et al. 2006b), or an apparent luminosity of L(>250 



GeV)= 2.6 x 10 jy erg s _l (adopting a distance of 3.8 Mpc). 

No significant variability has been found on time-scales of 
28 min, nights and months (moon periods). From the error on 
the nightly flux (~ 5% Crab), we estimate the sensitivity of 
the dataset for flares with a duration of one day. For a ~ 4 cr 
detection in the Cen A dataset a flux increase during a single 
night by a factor w 20 would be needed (a factor * 15 for 
~ 3 cr). This can be compared to the VHE flux variation of 
factor ~ 5 - 10 detected from M 87 on time-scales of days. 

The results have been cross-checked with independent anal- 
ysis and calibration chains and good agreement was found. 

3. DISCUSSION 

Figure [4] shows the spectral energy distribution of Cen A 
ranging from X-rays to the VHE regime. The flux mea- 
sured by H.E.S.S. is clearly below all previous upper limits 
in the VHE regime. Extrapolating the spectrum measured 
with EGRET in the GeV regime to VHE energies roughly 
matches the H.E.S.S. spectrum, though the softer end of the 
error range on the EGRET spectral index is preferred. The 
recently launched Fermi observatory should provide a more 
accurate spectrum in the MeV-GeV range soon. 

Several authors have predicted VHE emission from Cen A, 

36 To derive the energy spectrum, a looser cut on the distance to the source 
is used (6 1 < 0.03 deg 2 ) to increase the number of photons (the standard cut 
is6> 2 < 0.015 deg 2 ). 
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and more generally discussed VHE emission from radio 
galaxies. A first class of models proposed the immediate 
vicinity of the supermassive black hole as the region of VHE 
emission. Neronov & Aharonian | ( |2007| l and Rieger & Aha- 
|ronian (2008 ) proposed a pulsar-type particle acceleration in 
the magnetosphere of the sub-Eddington accreting supermas- 
sive black hole, which has been applied successfully to the 
low-luminosity radio galaxy M 87. In this m odel the relation 
Lie * ^bh is expected (see Eq. (12) in Rie ger & Aharonian] 
2008). Since the mass of the central black hole and the VHE 
luminosity from Cen A are ~ 10 times lower than the ones in 
M 87, such a model could possibly be applied to the current 
data on Cen A. 

It has also been proposed that, similar to the mechanism at 
work in other VHE blazars, VHE emission could be produced 
in the inner jet regions in radio galaxies ( |Bai & Lee]|2001| 
Chiab erge et al.|2001) . Given the jet angle of 15-80° towards 



for these protons, resulting in VHE y-ray emission. 



Kachelriess et al. (2008 ) considered possible UHECR emis- 



the observer (see e.g. Horiuchi et al. 2006 and references 
therein), the relativistic boosting in standard scenarios would 
be small. Ghisellini et al. ( 2005 ) proposed a two-flow type 
model (Sol et al.|1989[ l, with a fast spine and a slower, mildly 
relativistic sheath propagating within the jet, which has been 
successfully applied to M 87 ( |Tavecchio & Ghisellini | 2008 ). 
Their model for the SED of Cen A (see Fig. 3 in Ghisellini 
|et al.||20"05] ), which was mainly constrained by the available 
CGRO data, cuts off at lower frequencies than VHE, but could 
possibly account for H.E.S.S. data by refining the parameters 
(see also in this context Marc owTth et al.|1998| l. |Lenain et al. 
(2008) modeled the VHE emission of Cen A with a multi- 
blob SSC model. In this model, VHE emission is expected 
to take place in the broadened jet formation zone, where even 
for a large jet angle, a few emission zones can move directly 
towards the observer and Doppler boost the emission. Their 
model prediction for the VHE emission from Cen A is well 
compatible with the H.E.S.S. data reported here (see Fig. 7 in 
|Lenain et aLpOOgj 



More extended VHE emission may also be expected from 
Cen A. In this context, Stawarz et al. (2006 ) proposed that y- 
rays emitted in the immediate vicinity of the active nucleus are 
partly absorbed by the starlight radiation in the host galaxy. 
The created e ± pairs are quickly isotropized and radiate VHE 
y-rays by inverse Compton scattering the starlight radiation. 
The small size of the resulting isotropic pair halo (~ 4 arcmin 
in diameter) is fully consistent with a point-like source for 
H.E.S.S., but could be resolved by the future CTA (Cherenkov 
Telescope Array) 37 observatory. Stawarz et al.'s prediction for 
the VHE emission of Cen A within this model, resulting in a 
steady flux of a few 10~ I3 erg cirT 2 s and a photon index of 
~ 2.6 at the TeV energ y range, agrees well with the H.E.S.S. 
spectrum (see Fig. 6 in |Stawarz et al.|20 06 ). 

Furthermore, hadronic models have been invoked to predict 
VHE emission from radio galaxies. In this context, Reimer 
et al. ( |2004| ) proposed a synchrotron-proton blazar model for 



M 87 where the high energy component in the SED is in- 
terpreted in terms of synchrotron emission from either pri- 
mary protons or secondary and created in the inner 
jet. At larger distances and in the peculiar case of Cen A, 
ultrarelativistic protons with E p < 10 18 eV could be released 
within the kpc-scale jet, and would be then confined in the 
host galaxy by interactions with the interstellar magnetic field. 
Both the dust lane of Cen A and thermal gas of the merger 
remnant could act as efficient, dense targets of cold material 

37 http://www.cta-observatory.org/ 



sion from Cen A in the view of the CGRO observations and 
current limits provided by imaging atmospheric Cherenkov 
telescopes. Depending on the choice for the UHECR injection 
function in their model, some solutions proposed for the pre- 
dicted VHE y-ray spectrum are compatible with the H.E.S.S. 
results reported here (see Fig. 1 in Kachelri ess et al.|2 008). 

It has also been argued that VHE emission could originate 
from the outer giant radio lobes (e.g. Hardcastle et al.|2008[ ) 
lying at ~ 4° from the core. However, the VHE y-ray ex- 
cess presented here only matches the position of the core, the 
pc/kpc inner jets and the inner radio lobes. Only upper limits 
at VHE have been derived for the southern outer lobe ( |Kabuki| 
|et al.|2007) . 



Recently, Croston et al. (2009) reported the detection of 
non-thermal X-ray synchrotron emission from the shock of 
the southwest inner radio lobe. The edge of this lobe is lo- 
cated ~ 5' fro m the nucleus an d reveals edge-brightened X- 
ray emission ( |Kraft et al.|2003[ ). While the position is ~ 3cr 
away from the best fit position of the VHE excess, it is well 
within the upper limit of the extension. Studying the spatial 
variation of the spectral index in X-rays across this shock re- 
gion, Croston et al.| ([2009 ) constrained the high energy cut-off 
in the electron distribution to be y max ~ 10 8 . They inves- 
tigated inverse Compton scattering the starlight radiation and 
the CMB from high energy particles in this lobe and predicted 
a VHE emission well compatible with the H.E.S.S. measure- 
ment reported here (see their Fig. 8) for both seed radiation 
fields. This study would suggest Cen A is analogous to a gi- 
gantic supernova remnant (SNR). 

Besides the components of the AGN hosted in Cen A, one 
might wonder whether other sources in the host galaxy could 
be responsible for the VHE emission. For example, |Kraft 



et al. ( |2001[ ) detected more than 200 X-ray point sources in 
Cen A, and IACTs have detected about 60 VHE sources in our 
Galaxy, many of which are associated with SNRs and pulsar 
wind nebulae (PWNe). However, source ensembles, such as 
the sum of the contributions from SNRs/PWNe, would require 
an unrealistically large number of sources (assuming a typical 
luminosity of ~ 10 34 ~ 35 erg s _1 above 250 GeV per source). 

Further information on the VHE excess position, more de- 
tailed spectral shape and variability studies are required to dif- 
ferentiate between the different models. However, if the VHE 
emission is due to a misaligned blazar-like process - such as 
leptonic or hadronic emission from the jet - the proximity of 
Cen A makes it a very good laboratory to further investigate 
emission processes and jet physics in blazars. If the VHE 
emission originates from another process - e.g. UHECR in- 
teracting with the interstellar medium, or a SNR type process 
at the shock - this would be very exciting, giving new insights 
into the physics in the VHE domain. 

Cen A represents a rich potential for future VHE experi- 
ments. Our current data are at the edge of differentiating the 
possible emitting regions. With higher sensitivity (factor 10), 
better astrometric accurac y and angular resolu tion (e.g. ~ 5" 
and ~ 1', respectively) ( |Hermann~ et al. 2007), CTA would 
allow the localization of the site of the VHE emission, and, 
possibly, reveal multiple VHE emitting sources within Cen A. 
More generally, the detection of VHE emission from Cen A 
together with the detection of M 87 and the galactic center 
poses the question of whether VHE emission (y-ray bright- 
ness) might be a general feature of AGN. While the sensitiv- 
ity of current generation experiments is probably too low to 
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answer this question, one can hope that the CTA experiment 
will be able to detect a large enough sample of sources to shed 
some light on this issue. 
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